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Having thus described our invention, what we claim as new and desire to secure by 
Letters Patent is as follows: 



1. A method for digital printing represented as a matrix I, comprising: 

forming a sequence of matrices P^ with entries 0 or 1 where a 1 at some entry in some P^ 
represents that this pixel will be printed at stage k; 

constructing as sequences of matrices Ik with entries in [0, 1], so that Io= I; 

determining, when considering all pixels in Ik for all succesive values of k, a next pixel 
having a largest weight indicating that said next pixel is to be printed next , so that P^., differs 
from Pk by a zero at said next pixel in Pk being replaced by a 1 at the same position in Pk.,; 

printing said pixel; 

replacing the value of said pixel in Ik by a 0 thus forming Ik+i; and 

incrementing the value of k, until enough pixels have been printed to represent the overall 
darkness of I by the prmted image. 

2. The method according to claim 1, wherein with each printing of a pixel, an order of printing of 
remaining pixels is redefined. 
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3. The method according to claitn 1, further comprising reordering the pixels to be printed with 
each printing iteration. 

4. The method according to claim 1, wherein an original image is denoted as said matrix I = {I(i, 
j)}, where 1 < i <Himage and 1 < j <Vimage, and 

wherein each element I(i, j) of I is a real number between 0 and 1 where 0 represents 
"white", 1 represents "black" and intermediate values represent levels of grey. 

5. The method according to claim 4, wherein said matrix forming includes: 

setting k = 0 and defining Po as the constant Himage by Vimage zero matrix. 

6. The method according to claim 5, further comprising: 

with each iteration of redetermining, replacing k by k+1 ; and 
setting Ik = Ik.i and Pk = Pk-i except for a single pixel. 

7. The method according to claim 6, further comprising: 

determining whether there has been a change in V(i, j, Mk); 

computing all weighted averages W(i, j, Ik) not computed in earlier iterations and setting 
W(i, j, Ik) = W(i, j, Ik.,) if W(i, j, Ik) is not computed again, such that when k=l, all said weighted 
averages are computed, and for k>l, only said weighted averages which have been modified at a 
previous iteration are computed. 
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8. The method according to claim 7, further comprising: 

ordering all W(i, j, Ik) by decreasing order; and 

considering Max(k) as a larger value of all W(i, Ik)s. 

9. The method according to claim 8, further comprising: 

determining whether j, Ik) = Max(k). 

10. The method according to claim 9, further comprising: 

for all pairs (i, j) such that W(i J, Ik) = Max(k), replacing PkCi, j) = 0 with Pk(i, j) = 1 in 

Pk, and Ik(i, j) = I(i, j) with Ik(i, j) = I(i, j)-l in Ik. 

11. The method according to claim 9, further comprising: 

computing GreyTotal(k) = Zij Pk(i, j). 

12. The method according to claim 1 1, wherein if GreyTotal(k) < GreyTotal, then a next 
iteration is begun. 

13. The method according to claim 1 1 , further comprising: 

if GreyTotal < GreyTotal(k), and |GreyTotal-GreyTotal(k)| < |GreyTotal-GreyTotal(k-l)|, 

then set P=Pk. 
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14. The method according to claim 1 1, further comprising: 

if GreyTotal < GreyTotal(k), and |GreyTotal-GreyTotal(k)| > |GreyTotal.GreyTotal(k-l)|, 



then setP = Pk.i, 



15. The method according to claim 1, wherein a total number of black dots being printed is 
variable. 



1 6. The method according to claim 15, further comprising: 

providing a multi-tone printer, such that a lightest grey darker than Max(k) is printed, and 
Ik(i. j) = I(i, j) - 1 is replaced by 1^(1, j) = I(i. j) which represents a multibit grey level being 
printed. 

17. The method according to claim 9, further comprising: 

compensating for isolated pixels by making weights in V(i, j, M) depend on what is 
printed at (i, j). 

•1 8. The method according to claim 1 , further comprising: 

encoding the set of pixels to be printed after half-toning, in the matrix P={P(i, j)}, where 
U i ^ Himage and 1 ^ j ^Vimage, wherein each element P(iJ) of? has a value of either 0 or 1, 
where "0" represents "white" and "1" represents "black". 
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1 9. The method according to claim 1 8, further comprising: 

given a Himage by Vimage matrix M with entries in [0,1], selecting a neighborhood V(i, 
j, M) for each (i, j), whose shape and size selectively depends or not on (i, j), and a set of weights 
associated to all pixels in V(i, j, M). 

20. The method according to claim 19, wherein W(i, j, M) represents a weighted average of the 
elements of M in V(i, j, M), and GreyTotal= Zy I(i, j). 

21. A method of printing, comprising: 

forming a matrix of pixels; 

determining an order of printing of said pixels, said determining including finding a 
weight of said pixels and printing a pixel having a highest weight; and 

reordering the remaining pixels and printing a pixel having the greatest weight of the 
remaining pixels imtil all pixels have been printed. 

22. The method according to claim 21, wherein the pixels are printed in turn based on the 
darkness of the local image being printed. 

23. A method for digital printing, comprising: 

forming a sequence of matrices Pk with entries 0 or 1 where a 1 at some entry in some Pk 
represents that this pixel will be printed at stage k; 
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constructing as sequences of matrices Ik with entries in [0, 1], so that lo = I; 
determining, for a pluraUty of pixels in Ik for all successive values of k, a next pixel 
having a largest weight indicating that said first pixel is to be printed next; 
printing said pixel; and 

determining for each pixel of the remaining ones of pixels of said plurality of pixels a 
printing order of said remaining pixels such that subsequent pixels of said remaining ones of 
pixels of said plurality of pixels having a largest weight among the remaining pixels, are 
subsequently printed. 

24. A system for printing, comprising: 

means for forming a matrix of pixels; 

means for determining an order of printing of said pixels, said determining including 
finding a weight of said pixels and printing a pixel having a highest weight; and 

means for reordering the remaining pixels and printing a pixel having the greatest weight 
of the remaining pixels until all pixels have been printed. 

25. The system according to claim 24, wherein with each printing of a pixel, an order of printing 
of remaining pixels is redefined. 

26. The system according to claim 24, wherein said reordering means reorders the pixels to be 
printed with each printing iteration. 
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27. The system according to claim 24, wherein an original image is denoted as said matrix I = 
{I(iJ)}, where 1 < i <Himage and 1 < j <Vimage, and 

wherein each element I(i, j) of I is a real number between 0 and 1 where 0 represents 
"white", 1 represents "black" and intermediate values represent levels of grey. 

28. The system according to claim 27, wherein said matrix forming means includes: 

means for setting k = 0 and defining Po as the constant Himage by Vimage zero matrix. 

29. The system according to claim 28, further comprising: 

with each iteration of redetermining, means for replacing k by k+1 ; and 
means for setting Ik = Ik-i and Pk= Pk-i except for a single pixel. 

30. The system according to claim 29, further comprising: 

means for determining whether there has been a change in V(i, j, Mk); and 
means for computing all weighted averages W(i, j, Ik) not computed in earlier iterations 
and setting W(i, j, Ik) = W(i, j, Ik-i) if W(i, j, Ik) is not computed again, such that when k=l , all 

said weighted averages are computed, and for k>l, only said weighted averages which have been 
modified at a previous iteration are computed. 



31 . The system according to claim 30, further comprising: 
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means for ordering all W(i, j, Ik) by decreasing order; and 
means for considering ]V[ax(k) as a larger value of all W(i, j, Ik)s. 

32. The system according to claim 31, further comprising: 

means for determining whether W(i, j, Ik) = Max(k). 

33. The system according to claim 32, further comprising: 

for all pairs (i, j) such that W(i, j, Ik) = Max(k), means for replacing Pk(i, j) = 0 with Pk(i, 
j) = 1 in Pk, and Ik(i, j) = I(i, j) with Ik(i, j) - I(i, j)-l in Ik. 

34. The system according to claim 32, further comprising: 

means for computing GreyTotal(k) = Zi j Pk(i, j). 

35. The system according to claim 34, wherein if GreyTotal(k) < GreyTotal, then a next iteration 
is begun. 

36. The system according to claim 34, further comprising: 

if GreyTotal < GreyTotal(k), and |GreyTotal-GreyTotal(k)| < |GreyTotal-GreyTotal(k-l)|, 
means for setting P=Pk. 
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37. The system according to claim 34, further comprising: 

if GreyTotal < GreyTotal(k), and |GreyTotal-GreyTotal(k)| > |GreyTotal-GreyTotal(k-l)|, 
means for setting P = Pk-i. 

38. The system according to claim 24, wherein a total number of black dots being printed is 
variable. 

39. The system according to claim 38, further comprising: 

a multi-tone printer for printing, such that a lightest grey darker than Max(k) is printed, 
and Ik(i, j) = I(i, j) -1 is replaced by Ik(i, j) ^ I(i, j) which represents a multibit grey level being 
printed. 

40. The system according to claim 32, further comprising: 

means for compensating for isolated pixels by making the weights in V(i, j, M) depend on 
what is printed at (i, j). 

41 . The system according to claim 24, further comprising: 

means for encoding the set of pixels to be printed after half-toning, in the matrix P={P(i, 
j)}, where 1 < i < Himage and 1 < j <Vimage, wherein each element P(ij) of P has a value of 
either 0 or 1, where "0" represents "white" and "1" represents "black". 



42. The system according to claim 41, further comprising: 
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given a Hiiriage by Vimage matrix M with entries in [0,1], means for selecting a 
neighborhood V(i, j, M) for each (i, j), whose shape and size selectively depends or not on (i, j), 
and a set of weights associated to all pixels in V(i, j, M). 

43. The system according to claim 42, wherein W(i, j, M) represents a weighted average of the 
elements of M in V(i, j, M), and GreyTotal= Zjj I(i, j). 

44. The system according to claim 24, wherein the pixels are printed in tum based on the 
darkness of the local image being printed. 

45. A signal-bearing medium tangibly embodying a program of machine-readable 
instructions executable by a digital processing apparatus to perform a method of printing, said 
method comprising: 

forming a matrix of pixels; 

determining an order of printing of said pixels, said determining including finding a 
weight of said pixels and printing a pixel having a highest weight; and 

reordering the remaining pixels and printing a pixel having the greatest weight of the 
remaining pixels until all pixels have been printed. 
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46. A signal-bearing medium tangibly embodying a program of machine-readable instructions 
executable by a digital processing apparatus to perform a method of printing, said method 
comprising: 

forming a sequence of matrices Pk with entries 0 or 1 where a 1 at some entry in some Pk 
means that this pixel will be printed at stage k; 

constructing as sequences of matrices Ik with entries in [0, 1], so that Iq = I; 

determining, when considering all pixels in Ik for all successive values of k, a next pixel 
having a largest weight indicating that said next pixel is to be printed next , so that Pk+i differs 
from Pk by a zero at said next pixel in Pk being replaced by a 1 at the same position in Pk+i; 

printing said pixel; 

replacing the value of said pixel in Ik by a 0 thus forming Ik+i; and 

incrementing the value of k, until enough pixels have been printed to represent the overall 
darkness of I by the printed image. 
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